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mm MM TTTTTTTTTT HH HH DDDDDDDD EEEEEEEEEE Xx XX PPPPPPPP 
MMMM MMMM TT HH HH DD bd EE XX XX PP PP 
MMMM MMMM TT HH HH DD DD EE XX XX PP PP 
“MMM MM 1T HH HH DD D EE XX XX PP PP 
“M oOMM OM TT HH HH DD DD EE XX XX PP PP 
mM “mM TT HHHHHHHHHH = DD DD EEEEEEEE XX PPPPPPPP 
2M mM TT HHHHHHHHHH = DD DD EEEEEEEE XX PPPPPPPP 
“My mi TT HH DD DD EE XX XX PP 
mm nM TT HH HH DD D EE XX XX PP 
“m mM TT HH HH DD DD EE XX XX PP or 
4M a TT HH HH DD DD EE XX XX PP ere 
uM i TT HH HH DDDDDDDD  EEEEEEEEEE XX XX PP sot 
“Mm oN TT HH HH DDDDDDDD  EEEEEEEEEE Xx XX PP ; 
| 
LL WII! SSSSSSSS 
LL eeeeg SSSSSSSS 
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; File: MTHDEXP.MAR Edit: RNH1014 


~TITLE MTHSDEXP 
~ IDENT /1-014/ 
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COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


0 
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THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
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SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITA! 
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; FACILITY: MATH LIBRARY 
; ABSTRACT: 


MTHSDEXP is a function which returns the double floating point 
le Be bs bi floating point argument. 
e call is standard call-by-reference. 


VERSION: 1 


HISTORY: 
; AUTHOR: 


Peter Yuo, 15-Oct-76: Version 0 
Mary Payne and Jud Leonard, 0O3-JAN=78: Version 1 


MODIFIED BY: 
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Detailed Current Edit History  6=SEP-1984 11:22:02 CMTHRTL.SRCJMTHDEXP.MAR; 1 a 3) 


-SBTTL HISTORY ; Detailed Current Edit History 
; Edit History for Version 1 of MTHSDEXP 


1-1 Split single and double precision routines into two parts; 
Used gore Segerene and faster algorithms provided by M. Payne. 

Fixed bug causing unexpected snteger overflow. JMT 24-Mar-78 

Change from rational approx. to Chebyshev_approx. for 

more accuracy with overhang. MHP 18-May-78 

Detect underflow properly. TNH 16-June-7 

Remove extra instruction. TNH 16-June-78 

Move .ENTRY mask to module header. TNH 14-Aug-78 

Update version number and copyr ght notice. JBS 16-NOV-78 

4 onge symbols MTH__FLOUNDMAT and MTH__FLOOVEMAT to 

MTHSK_FLOUNDMAT and MTHSK_FLOOVEMAT. " JBS 07-DEC-78 

Add ''_” to the PSECT directive. JBS 22-DEC-78 

Declare externals. SBL 1 Pips Ob 

Use only through R6. SBL 21-Sept-1979 

Included logic for large arguments (where EMOD Loses 

significance). RNH 23-JUN-81 

- Changed W* to G* on calls th MTH$$SIGNAL and MTHSSJACKET_TST 

RNH 09-Sept-1981 
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: INCLUDE FILES: 


ete 


Sete Sete 
o 
z 
oO 
n” 


»PSECT 


OWN STORAGE: 
CONSTANTS: 


; Table to be used for a+ thy > These constants here have been 
er program ROOT16 on PDP-10. 


; 2**(0/16) = 1.0 


; verified by M 


TABHI: 


16-SEP- 


DECLARATIONS 


MTHJACKET 


EXTERNAL SYMBOLS: 


GBL 
MTHSK_FLOUNDMAT 
MTHSK~FLOOVEMAT 
MTHS$SSIGNAL 
MTHSSJACKET_TST 
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wm 
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PSECT DECLARATIONS: 


_MTHSCODE 


none 
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; Declarative Part of Module 


-MAR 


R3, R4, RS, R6> : 


-ENTRY mask + int oy 


3; Extension for operand in EM 


ro 


f enable 
DD 


; define SF$ (stack frame) symbols 


PIC,SHR,LONG,EXE ,NOWRT 


3 program section for math routines 


; 2e0(1/16) 
; 2**(2/16) 


2** (3/16) 
2** (4/16) 
2**(5/16) 
2** (6/16) 
2**(7/16) 
2** (8/16) 
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21446,°0027110 
21510,*0020320 
21730,°0037542 
22136,°C106532 
22044 ,*0137341 
21457,°0015020 
21746,°0154770 
22167,*0046406 
22002 ,*0127400 
21613, 0003526 
21540,°0167132 
21274, *0062600 
21410,*0031300 
21673,*0040264 
21661,*0167167 


6 
6 
6 
6 
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Oo Oo 20 20 200 


Or 
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: Constants used in evaluation of 


“00 $3 0,°001 
*00 4,°01 
“00 
“004 


4 
120,*001 
461,°011 
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2**(9/16) 
2**(10/16) 
2**(11/16) 
2**(12/16) 
2e* (13/16) 
2**(14/16) 
2**(15/16) 


DECIMAL: 0. 
; DECIMAL: 0. 
; DECIMAL:-0. 
; DECIMAL: 0. 
; DECIMAL: 0. 
; DECIMAL:-0 
DECIMAL:-0. 
; DECIMAL:-0. 
; DECIMAL:-0. 
; DECIMAL:-0. 
; DECIMAL:-0. 
; DECIMAL:-0. 
; DECIMAL:-0. 
; DECIMAL: 0. 
; DECIMAL: 0. 
; DECIMAL: 0. 


3; DECIMAL: 
; DECIMAL: 


dO 

2252169616881804D-17 
2712242510500122n-17 
5861402647731367D-17 
12064576472234940-16 


-8930877995013540D-17 


2373071989573779D-17 
62572408308818800-17 
1340620676392399D-16 
7084371812598154D-17 
37683790651871620-17 
3048384309613603D-17 
1276624235300040D-17 
1845830375854930D-17 
5075495866202897D-17 
48228430606756190-17 


polynomials - small arguments 


0.2480427857745020D-04 
0. 1984369200268758D-03 
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0000 0000 
00000009 


AA3B 4288 
17FO 295C 


C800F7D1 72173E31 
73003C76 357927EF 
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ONS ; Declarative Part of Modu 6-SEP-1984 11:22:02 (CMTHRTL.S ROOM THDEXP.MAR; 1 
«WORD  *0035 +200 554 
«WORD *0120720,°016147 : DECIMAL: 0.1388888879690042D-02 
“WORD  *0036410,°01 431 
<WORD *0072406,°017177 ; DECIMAL: 0.8333333262370290D-02 
<WORD 0 401 3052 
~WORD *01 3 “01 : : DECIMAL:  0.4166666666667950D-01 
-WORD  *00 401 : 2 
;WORD *01 335 “0000 ; DECIMAL:  0.1666666666667437D-00 
WORD — me DECIMAL: 0.5000000000000000D-00 
* WORD "6040177, 40177777 ‘ 
"WORD *0177777.°0177777 ; DECIMAL: 0.1000000000000000D+01 
“WORD *0040200-0 
WORD 0.9 : DECIMA 0. 10900000000000000+01 

DX°LN1=<.-DXPTB1>/ : no. + ontties in table 


PRIPIPIPININYNININININININIPYIPINPIPININIPINIPINININIPIPIPINININININPINININININIDY 2 2 OO OO es 


WNINIAIWINIMIAIUAIDIPININIPINININININ 2 2 9 PH FH MQOOOOOOOCOCOCOCOVOOOOOOOOOM® 
NAVE WN OC OOD NAUE WN O OO NA UE WWI OO OD NAU EWN OO OODNOU EWN O0 


; Constants used in evaluation of 


DXPTAB: .WORD “00 $261 .20061157 
WORD “0021 $7,201 5016 
“WORD *00 31 7,°0164 33 
“WORD *014 837,201034 3 
“WORD 0027 41,2010¢211 
[WORD “00 40212012 73 
WORD 700306 $.20141 77 
“WORD *003321 20147273 
“WORD *0032435.°011253 
“WORD 20076728, 70074645 
“WORD 2003414 **0054106 
[WORD *0134045-*0011314 
[WORD *0035565.°0176757 
“WORD 20176026, °0026166 
[WORD *0037061,°0071027 
[WORD *0173721.°0147572 
“WORD 0,0 
[WORD 0.0 
DXPLN=<.-DXPTAB>/8 
D_ 1eL062. -£: 
ORD *0041270,*0125073 
“0024534 .*0013760 
D_LN2_OV. er ML: 
QUAD  *XC800F7D172173E31 
D_LN2 OV. rr LO: 
QUAD  *X73003C76357927EF 


polynomial - regular args 


DECIMAL: 
DECIMAL: 
DECIMAL: 
DECIMAL: 
DECIMAL: 
DECIMAL: 
DECIMAL: 
DECIMAL: 


DECIMAL: 


0.3077130709430240D-15 
0.5682419384166091D-13 
0.9181219559808114D-11 
0.1271587192556359D-08 
0.1467610032291993D-06 
0.1355080777949815D-04 
0.93838479280898720-03 
0.4332169878499658D-01 


0 
no. of entries in table 


; LOG2(E) * 16 


; Hi 45 bits of ln2/16 


Low bits of |n2/16 


Lee, 


-( 
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-SBTTL MTHSDEXP = Standard Double Precision Floating EXP 


+ 

FUNCTIONAL DESCRIPTION: 

EXP = Double precision floating point function 

Uses a Chebyshev approximation, with overhang on Last step. 


CALLING SEQUENCE: 
Exponential.wd.v = MTHSEXP(x.rd.r) 
INPUT PARAMETERS: 


00000004 


00000004 ; define ee multiplier 


LONG = 4 

x = 1 * LONG ; Contents of x is the argument 
IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 

VALUE: floating exponential of the argument 
IMPLICIT OUTPUTS: none 
SIDE EFFECTS: 
Signals: MTH$_FLOOVEMAT if X > 88.028 with reserved ererend in RO/R1 (copied 
to the signal mechanism vector CHFSL_MCH_RO/R1 by LIBSSIGNAL). Associated 
message is: "FLOATING OVERFLOW IN MATH LIBRARY’. Result is reserved operand 
-0.0 unless a user supplied (or any) error handler changes CHFS$L_MCH_RO/R1. 


) 
MTHS_FLOUNDMAT if X =< “89.646 and caller has hardware enable Set. 
Hy Ate is set to +0.0. Associated message is: ‘FLOATING UNDERFLOW IN MATH 


NOTE: This procedure disables floating point underflow, enable ce oe 
overflow, causes no floating overflow or other arithmetic traps, an 
preserves enables across the call. 


407C -ENTRY MTHSDEXP, ACMASK ; standard call-by-reference entry 
r qiseete DV (and FU), enable IV 
a 


MTHSFLAG_JACKET 3 g that this is a jacket procedure 
MOVAB G*MTHSSJACKET_HND, (FP) 


$5 


6D 00000000°GF 9E 


COOOOOCOCOCOCOCOCOCOCOOOCCO OVO OCOOOCOCO COCO OCCOCOCOCOCOOOOSVOOCOOOCOOCOOCOOOOOoOOo 


81 3; set handler address to jacket 
4 3; handler 
81 8 ; in case of an error in special JSB 
81 9 3 routine 
50 06 BC 70 0181 390 MOVD ax (AP), RO RO7RI = user's arg 
1 10 B5 91 BSBB MTHSDEXP_R6 ; RO/R1 = special EXP(RO/R1) 


04 0187 292 RET ; return = result in RO/R1 


Oe anbaabys Praeinlong ges taeptaeenent at, TESEECIRRS QY:ISEAZ YAKNE NRESR OLB me. P20 cy 
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-SBTTL MTHSDEXP_R6 = Special DEXP routine 
Special DEXP = used by the standard, and direct interfaces. 


CALLING SEQUENCE: 
pove anything needed in RO:R6 


on 


z2v 


ro V04-00 Page 8 
MTHRTL.SRCJMTHDEXP.MAR; 1 (5). 


; input in RO 


whether enabled or not. 


REGISTERS USED: 
RO/R1 = floating argument, then result 
R2/R3 = temp 
R5 = integer scratc 


h 
R6 - integer part of X * LG2(E) * 16 (161+J) 


MTHSDEXP_R6:: ; special DEXP routine 


MTHSDEXP-R7:: Release 1 name 
Bicw3 


52 50 8000 8F 


AB #*x8000, RO, R2 
_ he 3E00 BF = AS 
38 0580 ef B81 
60 Ff 


Preliminary test for over/under flow 
R2 = 1st word of iX: 

SUBW3 3€00, R2, R3 R3 = 4 + unbiased exponent 
CMP BO, 88 


a*x 
#*x5B0, R3 Compare iXi with 

ST to more tests if LSSU 
else, -4 < unbiased ex 


10 
1) 
18 
14 
13 
1H 
15 
0 
1 BLSSU SMT 
; p<8 
Z no exceptions in EMODD or APPROX 
4280 8F 52 B61 5 CMPW R2, #*xX4280 ; Check for loss of significance in 
§ ; EMOD ( IX! >= Qed 
20. «6119 i BLSS EVAL 3; No loss of significance 
9 
0 
4 
40 
o] 
23 
44 
4 


; parts of X*16/ln2 must be obtained in seperate steps. 


; IX! >= 2006, EMOD will Lose significance so the interger and fractional 
. MULD3 D_16L0G2_E, RO, R2 


52 50 BB AF 65 Got tojoger part of X*16/Ln2 in 


CVTDL R2, R +J) as a longword and in 
5 R6, R R2/R3 in D format 
54 52 BB AF 6 MULD D_LN2_ OV_16_HI, R2, R4 Get fraction part of X*16/\n2 = 
16/ineat KX = (14J)#ln2/16 J 


R4, R 
D_LN2 OV_16_L0, R2 


0. 54 o¢ 
52 88 AF in RO/RI. 
0 52 62 SUBD 
50 Al AF % 
08 1 BRB 


RZ, R . 
MULD 0_{6L0G2_E, RO ; 
APPROX 3 
EVAL: EMODD 0D_16L0G2_E, #X_273, RO, R6, RO 

3; get X*16*LG2(E) with 

5 satoner part in R6 (=161+J) 

; ion in RO/R1 


frac 


TWH HHH HMMMMMOOOOCVCVOVTVTVTVFVVTCOCCOCCoCWWWVIOVVT TWO OVWOOOTOIOIOIIOIIION ve 


SSP FF > SVN OUMIS } WPI OOD DW W W © FTF 09.00.00 0 60 60.09.09 09 G9 0 490 09.09 0 49 OO OP ODED ODOT aD # 


$s 
$ 
94 
32 
i 
3 
0 Vv eee 
0 JSB MTHSDEXP_R6 
¢ return with result in RO/R1 
4 ; Note: This routine is written to avoid causing any integer overflows, 
5 ; floating overflows, or floating underflows or divide by 0 conditions, 
3 
09 
4 
FFB CF «608)— 50s rt} APPROX: POLYD RO,#DXPLN-1,DXPTAB ; use Chebyshev series 
4 ; with last coefficient 0 


PDI Se 4 9 — 4 2 4 4 4 4 — 4 —» > + 2 2s 0b 2s as ss 2s —» 2s as 5 > as 5 a» 2s as 5 ss as ns as 2s as tn 2 2 xe 


50 3 so that Last ADDD has overhang 


MTHSDE XP 
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R6 = Special DEXP rout :02 (CMTHRTL.SR RCM MTHDEXP.MAR; 1 ( 


55 56 FFFFFFFO 8F BICL3 #16, R RS 7; RS = J 
FOF3 CFG MULD —s- TABHI R .R : else MUL by 25064716) 
4343 CF4 4 ADOD TABLOLRSJ, R ; add in LO of 2**(J/16) 
FDE7 CF4 5 ADDD TABHILRSJ, R ; and then HT of 2**(J/16) 
a. ; pe ote Leer 
: = en 
50 6046 33 MOVAQ (RO)CR6], RO 3 pits i to Exe position estes | 
: ee nen 
OO7F BF 3 61 CMP RO, #*X7F ; test tor euertundert lew om | 
D $6 BLEQ EXCEPT ; see what exception is if ino or = 0 | 
a7 20$: RSB ; otherwise return result in | 
65 SMTST: | 
1 $6 BLSS 208 z exception if exp+4 > 11 
2400 8F 5 6 CMP ne #*x2400 z eliminate underflow from APPROX1 
0 : BLSS 108 : bypass if E**ARG = 1 | 
| 
0 ;+ 
71 ; Use Chebyshev series for small arg 


FEC6 CF O08 50 74 POLYD RO,#DXPLN1-1,DxXPTB1 ;Use Chebyshev series 
75 ; last term is 1; this will 
7 ; give desired overhang. 
Ae RSB 3 answer is OK, return 
50 6(08 79 10$: MOVD s*#1, RO ; Ex*exX is i, store it 
80 RSB 3 and return 
H | 
r : Handlers for software detected over/underflow conditions follow 
50 6 208: TSTD = RO : if big ARG > 0 goto OVERFLOW | 
28 7 BGEQ OVER 
§3 : Underflow; if user has FU set, signal error. Always return 0.0 | 
91 UNDER: 
52 3 MOVPSL Re : ~ = user's or jacket routine’s PSL 
00000000 ' GF 4 9 CALLS #0, G*MTHSSJACKET_TST ; = TRUE if JSB from jacket routine 
o% 0 94 BLBC R 3; branch if user did JSB 
52 AD 95 MOVZWL SFS$W_SAVE_PSW(FP), R2  : get user PSL saved by CA LL 
50 96 10$: CLRQ RO ; RO = result. LIBSSIGNAL will save in 
9 > CHFSL_MCH_ ROR so ony handler cen fixup 
0D 52 «06 9 BBC #6, R2, 208 3; has uSer énabled floating underflow? 
9 PUSHL (SP) i yes, return PC from special routine | 
7E 00° 400 MOVZBL #MTHSK_FLOUNDMAT, -(SP) ; trap code for hardware floa underflow | 
? 1 3 except on So ae _FLOUNDMAT ( in it VAX-11 | 
3 exception 
00000000'GF 02 4 ¢ CALLS #2, G*MTHS$SIGNAL 3 signal ka PC) | 
? : 208: RSB : urn 
4 $ EXCEPT: 
56 4 TSTL R6 ; test sign of 1; if 1 <0 


K 4 
MTHSDEXP Double Precision Floating Exponential 16-SEP-1984 01:1 AX/VMS Macro Vv04-00 Page 1 
1-014 MTHSDEXP_R6 = Special DEXP rout ine o7$Eb= 1984 91:38:03 :02 CMTHRTL.SRCIMTHDEXP.MAR; 1 2 (3) 
BLSS UNDER i go to underflow handler 


; Signal floating overflow, return reserved operand, -0.0 


Cee eres 
<= 
wm 
Po) 


else process for overflow 
return PC from special routine 
et floating overflow 


PUSHL (SP) 

MOVZBL #MTHSK_FLOOVEMAT, -(SP) 

ASHQ #15, #7, R ~ Tihtg et +f. operand 
e 

gnal mechanism vector (CHFSL _MCH_RO/R1) 


pe fe be fixed up by any error 


stonat” (condition, 
re 


LS #2, G*MTHSSSIGNAL PC) 
urn = RO restoredd from CHFSL_MCH_RO/R1 


Q0000000"GF 02 fe CAL 
RSB 


; l 
x 
dB 19 
E D 
re off 
50.01 F 9 


RoPoTnonono SS SS SS OOD 
FW OOOVNOUFS UT" O00 vu 


5 ek eh ee ee ee ek ee el ek ed Da 
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Symbol table 6-S 1:2 £03 MTHRTL. SR CIMTHDEXP.MAR: 1 ( 


ACMASK s 
APPROX R 02 
DXPLN = 
DXPLNI = 

XPTAB R 
DxPTB1 5 
D_16 0G _f R 
D-LN2_OV~16_H! bh 9 
D_LN2_OV_16_L0 R 

AL we 8 
m% ia 
MTHSSJACKET _HND eererere x f 
MTIHSSJACKET_TST eeeeeeee XK 0 
MTHSSSIGNAL” eeeeeeee =X 8 00 
MTHSDEXP Q0001A8 RG 0 
MTHSDEXP_R 0000188 RG 0 
MTHSDEXP_R 0000188 a” 0 
MTHSK_FLOOVEMAT eeeneene 0 
MTHSK_-FLOUNDMAT teeeeere x 00 
OVER 00000268 R 02 
SFS$W_SAVE_ PSW = 0000 es 
SMTS s'est BR 02 
TABH 00000000 R 02 
TABLO 00000080 R 8s 
UNDER 00000243 R 0 
x = 00000004 
X_273 = 00000088 

$oewn wen een een neces + 
: ; Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
,»§ 2 «= 444-44 : 69 9-3 NOPIC USR CON ABS LCL oe NOEXE NORD wnat NOVEC BYTE 
SABSS 0000 0:) 3° NOPIC USR CON ABS LCL NOSHR' EXE R&D T NOVEC BYTE 
_MTHSCODE b090027p 63 83 | PIC USR CON REL LCL SHR EXE RD NOURT NOVEC LONG 
_Debonebnnaenenneneenanena} 
: ; Performance indicators ! 

Phase Page faults CPU Time Elapsed Time 
Initialization 1 00:00:00.09 0:00:00.56 
ponnene processing 133 09: 0: 0:70 0:00:09.99 

ass :00:02. 00:07. 
pyabo), table sort a? 88 0 0:04 a 
Symbol table output 0:00: 0:08 8:80; 83 
Psect synopsis output 0:00:00. 0:00:00.0 
Cross-reference output 0:00: 3-8? 3:82:90 -09 
Assembler run totals 39 0:00:03.94 0:00:16.16 


The working sgt Limit was 1050 pages. : 
8608 bytes (17 pages) of virtual memory were used to buffer the intermediate code. 
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MTHSDEXP ; Double Precision Floating Exponential 16-SEP-1984 91:46:17 AX/VMS Macro V04-00 
VAX-11 Macro Run Statistics 6-SEP=1984 11:22:02 CMTHRTL.SRCIMTHDEXP.MAR; 1 


There were 10 pages of symbol table space allocated to hold 59 non-local and 5 local symbols. 
484 source Lines were read in Pass 1, producing 13 object records in Pass 2. 
9 pages of virtual memory were used to define & macros. 


tewrmececeme rasan aanreoccranecoen} 


Macro Library statistics ! 


' 
+ ween meee mewn ence e nen eso me a 


Macro Library name Macros defined 
_$255$DUA28: (SYSLIBISTARLET.MLB;2 pee ry Yr 
88 GETS were required to define 4 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHDEXP/OBJ=O0BJ$:MTHDEXP MSRC$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) +MSRCS: 


I 
R 
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